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Abstract 

We study the nonresonant three-body decays of B + — > K + it + and 5^ — > Ds*'' i"C°7r + . We 

(*) 

find that these decays can provide the information on the timedike form factors of Dg K. We 

(*)— 4- 

also explicitly investigate Bd — > D s K* + decays by discriminating the nonresonant contributions 
with the unknown Dg wave functions being fixed by the measured mode of Bd — ► D s K + . 
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Three-body decays of B meson have recently been noticed in the experiments by the 
analyses of the Dalitz plots and invariant mass distributions Jj. The study of the three- 
body decays provides not only the method to extract the CP violating phase angles ,2], but 
also the way to understand or search the uncertain particle states, such as /o(400 — 1200), 
/o(980), a (980) 3[ and glueballs Q. Moreover, it also helps us to build up the QCD 

approach for the nonresonant three-body decays H. As known that the charmless three- 

ft 

body decays of B meson are dominated by the so-called quasi-two-body decays jy], it is not 
easy to discriminate the nonresonant states from resonant ones. However, this may not be 
the case in those final states with charmed mesons. 

It has been demonstrated by Belle |.7j that — > K°K + decays are actually dom- 
inated by the three-body modes because the K°K + system is confirmed to be an J p = 1~ 
state by the analysis of angular dependence. The production of the decays can be thought 
easily by the consequence of b — > cud, while K°K + is produced by the created ss pair and 
ud. Therefore, if annihilation topologies are neglected, the dominant topologies for the de- 
cays correspond to Bd — ► D"~ and a outgoing pair of K°K + . The formers are described by 
B(i — > form factors, calculated by some QCD approaches such as perturbative QCD 

(PQCD), quark model, QCD sum rules and light-cone QCD sum rules, while the latters 
respond the times-like form factors which can be fixed via the connection to electromag- 
netic form factors and fitting with experiments pf. Although time-like form factors of KK, 
denoted by p l °-^ KK i are not easy to be formulated in theory, Bd — > D^*'~ K° K + is a good 
candidate to study the nonresonant three-body decays. According to the observations of 
Belle, we know that the branching ratios (BRs) of Bd — > D^*'~ K° K + decays are of O(10~ 4 ). 
. t means that the effects of F°^ KK are not small. Moreover, following the analysis of Ref . 
I|, we see that the peaks in the KK spectra of the decays locate at around 1.5 GeV. This 
can be understood that the region is actually governed by PQCD where the proper hard 
scale is around \/K.m B [5] with A = m B — m b . By taking A ps 0.48 and m B = 5.28 GeV, 
the value of Arris is 1.6 GeV, quite close to the consequence of Ref. 0. In some senses, 
PQCD approach can deal with the three-body decays by combing with the experimental 
fittings of time-like form factors. 

Inspired by the large BRs of Bd — > D^*'~ K° K + decays, one can speculate that the three- 
body decay related to form factors (DP\V^ — A^\B) can be also large, saying (9(10~ 5 — 10 -4 ), 
where D(P) is the charmed (pseudoscalar) meson and V^(A^) is the vector (axial-vector) 
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current. In theoretical viewpoints, the question is hard to answer since there are too much 
unknown form factors involved and no direct experimental data related to them. However, 
we still can give some conjectures on the relevant decays. Firstly, in terms of the concept 
of two-meson wave functions 10], the DP system could be described by a set of wave 
functions for (DP\Vfj, — A^\B) and they can be related to the time-like form factors of DP, 
denoted by F°^ DP . Therefore, if the K meson is massless particle, the threshold invariant 
mass, expressed by uj, to generate the KK pair is about 0. However, unlike the case of the 
KK pair production, since charmed meson is a massive particle, to produce the DP pair 
it should start from mo- By assuming that the peak of the DP pair spectrum is around 
mo + yAme, we can expect that the BR associated with F°^ DP form factors should be 
smaller than that associated with F°^ KK because the dominant form factors have been 
shifted to a larger to region and their values are small, compared to F°^ KK at yAm^. 
Moreover, if the third particle of the involving three-body decay is a light meson, although 
the allowed u of DP could reach the value of m# ( it is m# - mo in the Bd — > D^~ K°K + 
decays), due to the suppression of phase space factor (1 — u 2 /m 2 B ), the effects of the large ui 
are not important. Therefore, the available phase space is smaller than that in the decays of 
B d -»■ D^-R°K + . Hence we conjecture that the BR of B — > DPP' should be smaller than 
those of Bd — > D^~ K° K + decays, where DP system and B meson have the same light 
spectator. Note that the chosen examples have the same weak Wilson coefficients (WCs). 

Because there is no any direct information on F°^ DP ) in order to confirm our conjectures, 
we suggest that the observations of B + — > D^ K + tt + and Bd — > D^ K°tt + , illustrated 
by Fig. ^ can help us to find the answer. Since the former modes correspond to pure 
three-body decays and once they are measured in experiments, we immediately know what 
the effects of F°^ DP are. However, besides the nonresonant three-body decays, the decays 
of B d -> D { * ] ~K q k + also involve resonant states B d -> D { * ] 'K* + (K* + -> K°n + ). Recently, 
Belle ^lj and Babar have measured the relevant two-body decay Bd — > D~K + to 
be (4.61^ ± 1.3) x 1(T 5 and (3.2 ± 1.0 ± 1.0) x 10" 5 , respectively. One expects that the 

(*)— 4- 

Bd — > Ds K* + decays should have the same magnitudes in BR. Probably, the suggested 
three-body decays also have the same BRs in order of magnitudes. In order to understand 
more on the nonresonant parts in Bd — > D^ K°7r + , it is important to know how large the 
contributions are from quasi-two-body decays. In this paper, we want to make a detailed 
analysis on Bd — > Ds K* + decays. 
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FIG. 1: The topologies for the nonresonant three-body decays (a, b) B + — > D$*^ K + ir + and (c) 
Bd — > K°tt + . The dots denote the weak vertices. 

Since the considered decays correspond to the b — > cdu transition, we describe the effective 
Hamiltonian as 



q=u,c 



0[ q) = d aq/3 c p b a , O 



(1) 



where q a qp = q , Q 7 At (l — 7k )q/?, o>((3) are the color indices, V q = V* d V c b are the products of the 



CKM matrix elements and Ci^ip) are the WCs |l4j. Conventionally, the effective WCs 
of ci2 = C1 + C2/ 'N c and a\ = C2 + Ci/N c with N c = 3 being color number are more useful. It 
is known that the difficulty for studying exclusive hadron decays is from the calculations of 
matrix elements. In order to handle the hadronic effects, we employ the PQCD approach in 



which the transition matrix element is described 
and the hard amplitude of the valence quarks 
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i>y the convolution of hadron wave functions 
. Although the PQCD approach suffers 



singularities from end-point region, they could be smeared after the threshold and the fcy 
resummation effects are included. The latter arises from the introduction of the parton 



transverse momentum 
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181 ] - In the literature, the applications of PQCD to exclusive B 



B 



7] 
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meson decays, such as B — > B — > 7itt(KK) 
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2J, B -> <jm(K) 
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2J, B -> pK |25| decays and B -> K^U+l' have 



3een studied and found 
27, ba. We think that 



that all of them are consistent with the current experimental data 
the same approach can be also used to the considered cases here. 

It has been shown that by the reality of hierarchy mo >> m D (*) » A, the distribution 



amplitudes of B(Di*^) mesons can be described by |2£ 



(0|6((V(*),|B,pi) 



V2A, 



dxe 



-txpi -z 



c JO 



(2) 
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(D s ,p 2 \d(0) jC (z)i\0) = "^f / dxe ta »-'{'y B \fo + m Da ] li Q D .(x)} i 
(D:,p 2 \d(0)Az)i\U) = ^J^dxe^^Mh + rnn^n^x)}, (3) 



where e 2/1 is the polarization vector of D* meson, the normalizations of wave functions are 
taken to be dx& B , D (*-)Jx) = L^l'ij/V^c an d f B i D (*)\ are t ne corresponding decay 
constants. Although the decay constants and wave functions of D* between longitudinal 
and transverse polarizations are different generally, for simplicity, we assume that they are 
the same. Since the effects of transverse polarization parts are always related to the factor 
V2 = m B * s /M B , one expects that their contributions are much smaller than those from 
longitudinal parts. As to the distribution amplitude of the K* meson, we refer to the results 
derived from QCD sum rules 2iJ and summarize them as 

(0\u(0) jS (z)i\K*, P3 ) = -±= f ice^4m Jf .^] [j $ lr .(i) + M^(x) 

V ^ V c JO L 

+m K *[I]i j ® s K *(x) + m K *[hT]i3®K*{ x ) 

+[hTfa)ij® T K*(x) + rriK^^x)] , (4) 

where s 3 l{t) denote the longitudinal (transverse) polarization vectors of K* meson, C M = 
a _/p ?J ■ n_ in which n_ is parallel to the large component of p 3 , and & K l cor- 
respond to twist-2 wave functions while the remains stand for the twist-3 ones. In ad- 
dition, in terms of light-cone coordinate, the momenta of various mesons and the light 
valence quarks inside the corresponding mesons are assigned as: pi = m B /v / 2(l, 1, 0t), &i = 
m B /V2(x 1 ,0,k 1T ); p 2 = m B /V2(l, rf, T ), k 2 = m B /V2(x 2 ,0,k 2T )] P3 = m fl /\/2(0,l - 
r 2 2 , T ), k 3 = m B /V2(0, (1 - r 2 2 )x 3 , k 3T ). 

As usual, we have the decay rate for B c i — > D~K* + decay as 



r = — —^-\a, - 

16tt 



L D~K*+\ > 



where P c = \p 2z \ = \p&s\ an d the amplitude of A D - Rt+ is given by 

A D7K *+ = IbF d - k , + + M D7K . + . (5) 

The first (second) term in Eq. (j3J) comes from the factorizable (nonfactorizable) contribu- 
tions, illustrated by Fig. |2t (J2b)- The hard amplitudes F D - K , + and M D - Rt+ are expressed 
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FIG. 2: The topologies (a) factorizable (b) nonfactorizable effects for the decays Bd — > K* + . 
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DTK" 



by 

pi poo 

-8nC F m 2 B / d[x] / [b]d[b]$ Ds (x 2 ) 
Jo Jo 

| [ (r 2 _ (i _ 2r 2 > 3 )$^(C) + r 2 r 3 ((l - 2z 3 )$U0 + (1 + 2x 3 )$U0 
xEl(tl) + [(1 - r 2 2 )x 2 $^(C) - r 2 r 3 (l + x a )$^(C)]^(0}, 

lQnC F m 2 B ^2N c d[x) [b]d[b)^ B (x 1 ,b 1 )^ Ds (x 2 ) 
Jo Jo 

{ [ - (r 2 2 x 2 + (1 - 2r 2 )x 3 )$K*(C) + r 2 r 3 ((x 2 - x 3 )& K *(() 

+ (X 2 + * 8 )*5c*(C))] + [fo - HD<*MC) 

-r 2 r3((x3-x 2 )$^(C) + (2 + x 2 + X3)$^(C))]^ 2 (tJ)}, 

with r 3 = m^*jm B and £ = 1 — x 3 . In our considerations, the small effects from r\ 
and r 2 ~ 0.05 are neglected. The evolution factors £ l a rt\ are defined by 



Ci(0 + 



c 2 (4; 



a,(4)5. 



(f f K({*},{6}), 



f)OB+D 3 +K*^f) n f 
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where £ denote the hard scales and are chosen as 



1 1 



1 1 



*i = max{m B ^J (1 - rf> 3 , — , — ), t\ = max{m B ^J (1 - rf)x 2 , — , — ) 

t J f = max(m Bx fFf, m B J (1 - r%)x 2 x 3 , — , — ), 
F 2 = (1 - r|)(x x - x 2 )x 3 , Fl = xi + x 2 + (1 - rf)(l - ai - x 2 )x 3 . 



(6) 



(7) 
0.028 



(8) 



Here, SW*+d s = Sjf.^ ,p 3 -k 3 )S Ds (k^) and S B+Ds+K * = S B (kf)S Ds (k^)S K *(k 3 ,p 3 -k 3 ) 
are the associated Sudakov factors. We note that only the light valence quarks of B and 
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mesons have the Sudakov effects :lf|. h a j describe the hard functions arising from the 
propagators of gluon and internal valence quark. Their detailed expressions with threshold 
resummation effects can be found in Ref . |18[ . 

It is known that besides the longitudinal polarization, there also involve two transverse 
polarizations in B — > VV decays [31], where V denotes the vector meson. Therefore, the 
Bd — > D*~K* + decay amplitude will be more complicated than those in B — > VP or PP 
decays with P being the pseudoscalar meson. In terms of helicity basis, the B& — > D*~K* + 
decay rate is written as 

r = G2pPc V M m M {h \ 

B h=L,T 

where the superscript h denotes the helicity states of the two vector mesons. The amplitude 
M.( h > is decomposed into 



KA{h) _ * * 



a gn» + b . tf p » + i 1 e ^p 2aPw 

m D *m K * m D *m K * 



= m 2 B M L + m 2 B M N e* 2T ■ e* 3T + iM T e a ^ p e* 2Ta e; Tp p 2 ^p 3p , (9) 

with the convention of tr(7 5 /t J) j6 fi) — — 4ze a/37 ^ ^ajjpc^dp and the definitions of 

b 

m\M L = ae* 2L - e* L + e* 2L ■ pi e* 3L ■ p u 

m D *m K * 

m 2 B M. N = a e* 2T ■ e| T , 

M T = • (10) 

m Dt m K * 

Similar to Bd — > D S K* + , the decay amplitudes corresponding to each polarizations can be 
written as 

Ml(N,T) = fBFL(N,T) + M L ( NiT ), 

where the first (second) term expresses factorized (nonfactorized) effects. Each hard ampli- 
tudes can be formulated by 

pi poo 

F L = -8nC F m 2 B / d[x] / bd[b]$ D *(x 2 ) 
Jo Jo 

x{ M + (1 - r 2 2 )z 3 )<M0 + r 2 r 3 ($^(0 - $^(C))] 

- [(1 - r 2 2 )x 2 $ K *(0 + 2r 2 r 3 (l - x 2 )$U0] Oa)}, (H) 
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p 1 roo 

-8nC F m 2 B r 2 r 3 / d[x] / bd[b]§ T D »(x 2 ) 
Jo Jo 

X { (1 + X3 )*h (C) + (~^*K* (0 + (1 - *3«. (C)) 

- [(i + ^UO - (i - ^^(OKX)}, 



/*1 /*Q0 

-167rCVm|r 2 r 3 / / 
Jo Jo 



X 



;i - ^3)^(0 + ( + (i + x s )^(o 



mi) 



(12) 



(13) 



Mr 



167iC F MlJ2N ( 



x 



+ 



d[x] / [6]d[6]$ B (a:i, 61)^(^2) 
Jo 

1 - 2r 2 )x 3 $ x *(C) + r 2 r 3 ((x 3 + x 2 )$UO + (x 2 - x 3 ».(C)) f ,'(/ 
( r 2 + (1 _ 2r 2 2 )x 2 )$^(C) + r 2 r 3 ((2 - x 2 - x 3 )^»(0 



+ (x 2 -x 3 )$^(C))]^ 2 (tJ)} 



(14) 



/*1 /*oo 

M,v = -327rC F m|v / 2iVcr 2 r 3 / / [b]d[b)$ B (x u b 1 )® D *{x 2 )$ v K *{C)£'j{t)), (15) 

Jo «/ 

M T = -Q4nC F m 2 B y/2N c r 2 r 3 / / [b)d[b]^ B (x 1 , b 1 )^ D ,{x 2 )^ a K ,{C)£'j{t)). (16) 

Jo Jo 

From the above equations, we see clearly that the hard amplitudes for the transverse polar- 
izations are proportional to factors of r 2 r 3 and r|. One of two r 2 in the latter comes from 
the mass of charm-quark m c that we already set m c « m D w due to m^c*) — m c ~ A. 
To get numerical estimations, we model the S-meson wave function to be 



$b(x, &) = iV£X 2 (l - x) 2 exp 



[-5 




) 2 









(17) 



where b is the conjugate variable of the transverse momentum of the light quark, N B can 

be determined by the normalization of the wave function at b = and u B is the shape 

parameter. Since the wave functions of the K* meson are derived in the framework of QCD 

I 1 

sum rules, we express them up to twist-3 directly by [30] 



v/2A7 



:X{1 — X) 



1 + 0.57(1 - 2x) + 0.07C 2 3/2 (1 - 2x) 
8 



& K *(x) = -^=(o.3(l - 2x) [3(1 - 2xf + 10(1 - 2x) - 1] + 1.68C 4 1/2 (1 - 2x) 
+0.06(1 - 2x) 2 [5(1 - 2xf - 3] + 0.36 [1 - 2(1 - 2a;) (1 + ln(l - x))\ J , 

$^(x) = 2 ^g_ |3(l - 2x) [1 + 0.2(1 - 2x) + 0.6(10a; 2 - 10a; + 1)] 
-0.12a;(l - x) + 0.36[1 - 6a: - 2 ln(l - a;)] 1 , 



3/ 



A'* 



x(l - a;) 1 + 0.6(1 - 2a;) + 0.04C 2 3/2 (1 - 2x) 



1 + (1 - 2a;) 2 + 0.44(1 - 2xf + 0.4C 2 1/2 (1 - 2x) 



$ a K .(x) 



2^/2N~ c {4: 

+0.88C 4 1/2 (1 - 2x) + 0.48[2a; + ln(l - x)] 

4 ^_ |3(1 - 2x) 1 + 0.19(1 - 2x) + 0.81(10x 2 - 10.x + 1) 
-1.14x(l -x) + 0.48[1 - 6a; - 2 ln(l - a;)]} , 

with the Gegenbauer polynomials, 

Cl'\i) = \{Ze - 1) , C 4 1/2 (0 = i(35^ 4 - 30£ 2 + 3) , Cf{i) = ^ 2 - 1) . 

Although there are theoretical errors on the coefficients of the wave functions, we find that 
the allowed errors in these derived wave functions change the BRs only at a few percent 
level so that we will not discuss them further. In the PQCD approach, the wave functions 
represent the nonperturbative QCD effects and have the universal property. In principle, 
they can be determined by some specific measured modes. Hence, the unknown ub can 
be fixed by such as B — > (f>K decays, in which and K meson wave functions have been 
determined in the literature. Consequently, the remaining uncertain wave functions are the 

r»(*) 

D s mesons. 

As known that the BR of Bd — > D~K + has been observed by Belle and Babar so that we 
can take it as the criterion to model the relevant Dj wave functions; and because the mass 
difference between D~ and D*~ is not much, for simplicity, except the decay constants, their 
wave functions are taken to have the same behavior. It is worthful to mention that as the 
6-dependence on the wave function of B meson, for controlling the size of charmed mesons, 



we also introduce the intrinsic 6-dependence on those of charmed mesons. Hence, we use 
the wave functions of as 



4« 



6x(l — x)[l + cd 3 (1 — 2x)] > exp 



u Ds o 



(18) 



2y/2N c 

where cjj a and ud s are the unknown parameters. Although cd s is a free parameter, it can be 
chosen such that the meson wave functions have the maximum at x ~ A/m r M ~ 0.3 



for m c = 1.4 GeV. And then, we can fix u;d s through the observation of the Bd — > D~K + 
decay. 

In our numerical calculations, the values of theoretical inputs are set as: uog = 0.4, fg = 
0.19, / ( .) = 0.22(0.24), f$ = 0.20(0.16), f K = 0.16, m B = 5.28, m w = 1.969(2.112), 



m^* = 0.892 and = 1.7 GeV. With these values, we get F B ^ K = 0.35. Since the 
decay B d — > D~K + has been studied by Ref. j^, we use the derived formulas in Ref. 
directly. Although the decay Bd — > D*~K + has not been considered yet, due to the results 
only related to the longitudinal part of D*~ which should be similar to the longitudinal 
contributions of the D*~K* + mode, we also include it in our discussion. In Table HI we 
present the magnitudes of hard amplitudes. From the table, we see that the nonfactorizable 



TABLE I: The hard amplitudes of B d -> K* + (in units of 1CT 3 ) with c Ds = 0.9 and 

uo Da = 0.45. 





M D~K* + 


f B F L 


M L 


-0.088 -i 0.125 


2.37 + i 8.29 


-0.42 - i 0.09 


0.51 - i 5.12 


f B F N 


M N 


f B F T 


M T 


-0.01 + i 0.15 


0.22 -i 0.35 


0.23 -i 0.16 


-0.10 - i 0.24 



contributions of longitudinal parts are larger than factorizable ones. This could be under- 
stood by the mechanism of chirality suppression. The more remarkable example happens 
in the decays B — » titt. Because two 7r mesons are identical particles, under 22 ^ £3 trans- 
formation, the decay amplitude should be the same. However, from the same topology of 
Fig. |2t, we know that the four-momentum of the internal quark in one pion is opposite in 
sign to that in the another pion. That means the contributions from hard gluon exchange 
in opposite meson sides are canceled each other. Although Di*^~ and K* + are not identical 
particles, the cancellation should still exist. Since the momentum carried by 6-quark in 
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nonfactorizable parts, illustrated by Fig. Eb, is Pi — &i — ^2 — ^3 arid the light quark is 
k\ — hi — ks, the similar cancellations in nonfactorizable effects are not significant. From 
the Table EJ we also see that the transverse effects of B d — > D*~K* + are much smaller than 
those from longitudinal contributions. 

As a result, the predicted BRs of B d — > decays are displayed in Table HP (jlll|) 

for various values of cd s (wd s ) and ujd s = 0.45 (cd s = 0.9). To be more clear, we also 
show the behavior of the D s wave function with different values of cd s in Fig. El From 
the figure, we see that the larger cn is, the closer maximum of $ n (*) is to x = 0.3. Since 
the maximum of the D s wave function is located around 0.3 as shown early, a large cd s 
is preferred. Consequently, the BR of B d — > D~K + intends to be also large. However, it 
is necessary to introduce the 6-dependence in the wave functions to make the BR be 
more adjustable because the accuracy of experimental data is not good enough. According 
to our predictions, we find that the BRs of B d — > D*~K^ + decays are smaller than those of 
D~K^ + modes. Moreover, if the theoretical inputs are taken to fit with the observed value 
of B d -> D-K+, the BR products of BR(B d -> D ( S * ] K* + ) x BR(K* + -> K°ir+) should be of 
O(10" 5 ), with Br(K* + -> K°tt + ) = 2/3 Br(K* + -> (i^7r) + ). Therefore, if the observations 
7r + are close to O(10 4 ), the corresponding nonresonant three-body decays 
will be dominant. However, if experimental data conclude that the BRs of considered three- 
body final states are larger than our predictions significantly, but not approaching to 10~ 4 , 
it means that the pure three-body and quasi-two-body decays are comparable. On the 
other hand, if our predictions are consistent with the measurements of experiment, one can 
conclude that the three-body decays via the topologies of Fig. [T] are not important and 
negligible. 





TABLE II: 


The BRs (in 


units of 10 5 ) with the 


various 


values of cd s 


and ojd s = 


0.45. 




B d - 


- D-K+ 


B d -» D*-K+ 


B d 


-» D~K*+ 


B d - 


♦ D*~K* + 


0.9 




2.77 


1.34 




4.34 




2.00 


0.7 




2.33 


1.09 




3.79 




1.66 


0.5 




1.94 


0.87 




3.32 




1.39 



Finally, some remarks are given. In this paper, we assume that the mechanism of 

n 

the color transparency |3_3j] still dominates in B decaying to charmed mesons so that we 
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TABLE III: The BRs (in units of 10 5 ) with different values of ujd s and cd s = 0.9. 



uo Da B d —* DjK + B d ^D*-K+ B d ^DjK*+ B d — >■ D*~K* 



0.5 


2.34 


1.19 


3.79 


1.74 


0.4 


3.19 


1.50 


4.85 


2.30 


0.3 


4.18 


1.82 


5.98 


2.88 




FIG. 3: The wave function of <$>d 3 with cd„ = 0.9 (solid line), cd s = 0.7 (dashed line), and 
cd 3 = 0.5 (dashed-dotted line). 

don't have to consider the rescattering effects. From our results, we know that the PQCD 
approach can match the experimental data of the B d — > D~K + decay. With the same 
approach, our predictions on other decays, arising also from annihilation topologies, should 
be reliable. We note that the relevant charged B decays of B + — > Di*^ + K^ are also 
governed by annihilation contributions. Due to the suppression of CKM matrix elements, 
V u h V* d , the BRs are estimated to be of 0(1CT 8 ) [32]. Since the corresponding nonresonant 
three-body final states Di*^ + K~n + are also suppressed by the same CKM matrix elements, 
illustrated by Fig. 0J B + — > K~ ix + decays cannot give us more interesting information 
on three-body decays. 
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FIG. 4: Topology for nonresonant three-body decays B + — > Dg K 7r + . 
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